Proving the Generation of Genuine Multipartite Entanglement in a Single-Neutron 

Interferometer Experiment 
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We present experimental evidence of the generation of distinct types of genuine multipartite entan- 
glement between three degrees of freedom (spin, energy, and path) within single-neutron quantum 
systems. This is achieved via the development of new spin manipulation apparatuses for neutron 
interferometry and the entanglement is detected via appropriately designed and optimized non- 
linear witnesses. Via the applied criteria we reveal even finer properties of the type of the genuine 
multipartite entanglement that is produced in the experiment. In a subsequent analysis we show 
the extraordinarily high fidelity of the generated entangled states, proving the outstanding level of 
controllability of these systems. 

PACS numbers: 03.75.Dg, 03.67.Mn, 42.50.Dv, 03.65.Yz 



INTRODUCTION 

Entanglement is at the heart of the quantum infor- 
mation and communication technology It has been 
found in many distinct physical systems, e.g., photons, 
atoms, ions, neutrons, and even in quark- antiquark sys- 
tems of the second and third generation (e.g. [2}). The 
simplest entanglement is achieved with a bipartite sys- 
tem, while multipartite entanglement has been shown to 
enable different applications, e.g^ secure communications 
between many parties (e.g. [3j-|5[) or the basis for mea- 
surement based quantum computation (first introduced 
in Q). In addition, it may play a fundamental role in var- 
ious complex physical systems ranging from condensed 
matter systems (e.g. [7H10() to possibly even biological 
ones (e.g. Here, it is worth noting that entangle- 

ment of multipartite systems is essentially different from 
entanglement of bipartite systems. Multipartite entangle- 
ment, with which this paper deals, can e.g. be used for 
the study of quantum contextuality [l2| . whereas entan- 
glement between different particles (multiparticle entan- 
glement) can also be used to study quantum non-locality. 

In order to explore further sophisticated applications it 
is crucial to have a well founded understanding of multi- 
partite entanglement as well as to find systems in which a 
controllable physical implementation is feasible. A single- 
neutron system provides an extraordinarily controllable 
three qubit system and as such is perfectly suited to study 
the structure of tripartite Hilbert spaces and all predic- 
tions for correlations therein. We demonstrate multipar- 
tite entanglement by means of this system and in princi- 
ple open up a path to high precision tests of three qubit 
Hilbert spaces. For known physical applications that re- 
quire genuine multipartite entanglement, as e.g. quan- 
tum computing, intra-particle entanglement is probably 
not the first choice, however, for purposes of fundamental 



tests, which we aim for in this paper, intra-particle entan- 
glement serves just as well. In terms of potential physi- 
cal applications most progress has been made in photonic 
systems, as the techniques for manipulating photons lo- 
cally is very well developed 13]. Photons move at the 
speed of light, making them ideally suited for long range 
communication (such as quantum secret sharing). Nev- 
ertheless, although some attempts at photonic quantum 
computer have successfully been made, it is not trivial to 
hold and store photons at a place, which is an essential 
requirement for quantum computing of higher classes. In 
that respect other approaches have been pursued ranging 
from ion traps |14l| to proposals for quantum optomechan- 
ical systems [15| . 

Most physical implementations of multipartite entan- 
gled systems explore only one kind of internal degree of 
freedom that is entangled to the same internal degree 
of freedom of other physically disjoint systems at differ- 
ent locations. Another kind of entanglement is reported 
by entangling photons in all possible degrees of free- 
dom, which is called hyper-entanglement [16| . Entangle- 
ment purely between degrees of freedom in a single quan- 
tum system is extensively studied by using neutron opti- 
cal techniques in particular neutron interferometry |17| : 
starting with a demonstration of the violation of a Bell- 
like inequality [18], the implementation of a triply entan- 
gled Greenberger-Zeilinger-Horne (GHZ)-like state is re- 
ported [31 . Quite recently, a test of an alternative model 
of quantum mechanics a la Leggett [2(| as well as high- 
efficiency manipulations of a tripartite-entangled system 
[2ll ] are reported by using another strategy, namely, neu- 
tron polarimetry. It is worth noting here that these neu- 
tron optical techniques exploit matter- waves and are very 
suitable to study properties of multipartite entangled sys- 
tems of matter-waves: the efficiency of manipulating de- 
grees of freedom in a single neutron is rather high, which 
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enabled tests of alternative models of quantum mechanics 
with matter-waves. 



THEORY 

Multipartite entanglement is defined -as in the bipar- 
tite case- via an exclusion of separability. Generally, 
an n-partite pure quantum state is called fc-separable 
iff it can be written as a product of fc substates, i.e., 
I^fc-scp) = \<f>i) <8> 1^2} ® • • ■ <8> |0fc)- For mixed states there 
exists a straightforward definition, namely, a mixed state 
is fc-separable iff it has a decomposition into fc-separable 
pure states, i.e., p fc _ sep = EiPiK-sepX^LsepI- Entan- 
gled states consisting of more than two parties can be 
partially separable: it may be that they only are separa- 
ble with respect to certain partitions, while for other par- 
titions they are entangled. An important class of states 
pertains to genuinely multipartite entangled states, which 
are not biseparable (2-separable) with respect to all pos- 
sible partitions. A biseparable state is one that can be 
decomposed into a convex set of biseparable pure states. 
It is worth noting here that the partitions where the de- 
composed states are biseparable do not need to be the 
same for each decomposition: e.g., the tripartite state 
P = PPi <£> P23 + (1 - p)pi2 ® P3, with p G [0, 1], is bisepa- 
rable, though there is no specific bipartition of the whole 
state (with respect to which the state is separable). Nec- 
essary and sufficient criteria to detect biscparability (i.e., 
entanglement) are only known for the simplest bipartite 
case. Multipartite entanglement, on the other hand, is 
much richer in variety and more difficult to detect, since it 
does not need to be biseparable with respect to a specific 
decomposition. Various techniques to detect multipartite 



entanglement have been developed (e 



g. m 



22j-|25|). We 



adapt a specific technique from [23W25[. 

Two famous examples of genuine multipartite entan- 
glement, which have been shown to be very different in 
nature (e.g. Ref. [!§]), are the GHZ-state and the W- 
state (introduced in Ref. |26[ and 27[). These are the 
two states for which we present the experimental proof of 
production in the tripartite neutron system via multipar- 
tite detection criteria. In computational notation they 
read (a, b, c, d, e. . . amplitudes to be determined) 



\W) = a 1 101) +b e" /2 |011) +c 002} 



—\wi) 



~\w 2 ) 



= \w 3 ) 



\GHZ) = d e i7r/2 |101) + e|010) 



(1) 



(2) 



where the three subsystems denote three different prop- 
erties of single neutrons. The first subsystem denotes 
the path in the interferometer (path degree of freedom), 
which we also denote as path I (|I) = |0)) and II (|II) = 
1)), the second denotes the spin of the neutron (which 
can take the two values ±-| and we denote as | \) = |0) 



and I t) = and the third subsystem refers to the total 
energy (potential energy in an applied magnetic field plus 
kinetic energy) of the neutron (where the neutron can be 
at different total energy levels, denoted as \Eq) = |0), 
\E Q - hcu) = |1> and \E - 2hu) = |2}). 

For the detection of these two states we employ two 
distinct nonlinear entanglement witnesses, which are op- 
timally suited to detect GHZ- and W-like entanglement. 
The GHZ witness inequality employs a permutation op- 
erator n^, permuting the i th subspace of the two copy 
Hilbert space, reads 

Ighz ■= (3) 
3 

|(oio|p|ioi)| - ^(oioiofin^n.iofofof) < o, 

i=l 

and is satisfied for all biseparable states and maximally 
violated for the above defined GHZ-state. For the W- 
state a similar derivation leads to the witness inequality 



Iw '■= 



(4) 



^ \{Wi\p\Wj)\ - X! V ( w i w j\ U iP® 2Il i\ w i w 3) < °> 

which is again satisfied for all biseparable states and max- 
imally violated for the above defined W-state [with \wi) 
as defined in (JT}]. Thus any violation of any of these two 
inequalities is a clear proof of the presence of genuine mul- 
tipartite entanglement within the system. Furthermore 
both witnesses were recently shown to provide a good 
lower bound for an entro pic measure of genuine multi- 
partite entanglement (28l . l29j . To detect fc-separability, 
we employ the following inequality, which is satisfied by 
all fc-separable states [241 ]: 



h-seM ■ = VWlP^Pt^lW) 



(5) 



/ fc \ 2k 

{a} \»=1 / 

where |$) = \<j>i) (g> Ifo) is an arbitrary fully separable 
state, P ai is a permutation operator permuting the a^-th 
elements of |(?!>i) and \<j>2), Ptotal wholly permutes the \(j>i), 
and the sum runs over all fc-partitions {a}. 



EXPERIMENT 

Preparation of W- and 
Greenberger-Zeilinger-Horne-like states 

The experiment, illustrated in Fig. [TJ was carried out 
at the perfect-crystal neutron optics beam line S18 at the 
high flux reactor of the Institute Laue-Langevin (ILL) 
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[30[. By means of a Si perfect-crystal monochromator, 
a neutron beam with a mean wavelength of Ao = 1.92 A 
(AA/Ao » 0.02) was selected. Then the beam was po- 
larized vertically by magnetic-prism refractions. Because 
of the angular separation of the two sub-beams, only up- 
spin neutrons | j*) met the Bragg condition at the first 
interferometer plate. After this the beam, with an aper- 
ture of 4x8mm 2 , entered a skew-symmetric silicon perfect 
crystal neutron interferometer (IFM), which was adjusted 
to give a 220 reflection, thus coherently splitting the in- 
cident beam into the two spatially separated paths I and 
II at the first IFM plate. A parallel-sided Al plate was 
used as a phase shifter to vary the relative phase \ for 
the path degree of freedom prior to the coherent recom- 
bination of the two paths at the last IFM plate. A pair 
of water-cooled Helmholtz coils produced a fairly uni- 
form magnetic guide field BqZ of 2.3mT over the region 
of the polarized neutron beam. This field also defined a 
potential energy (Zeeman energy) of the neutrons, while 
the total energy of a neutron is equal to the sum of its 
potential energy and kinetic energy. Our energy degree 
of freedom is the total energy, as already utilized in our 



previous experiment [31] . To manipulate the neutrons, 
we used radio-frequency (RF) spin flippers, which change 
both the spin and the total energy, and a direct current 
(DC) spin flipper, which changes the spin without chang- 
ing the total energy. The change in total energy involved 
by a flip from spin-up to spin-down due to an RF coil 
operating at frequency w is given by — foui. By a tt flip 
we mean a flip at that most neutrons' spins are flipped 
(according to the efficiency of the spin flipping device), 
whereas for a ir/2 flip there is only a 50% probability for a 
neutron to be subject to a spin flip. Each time a neutron 
suffers a spin flip, it also suffers a phase shift e l7r / 2 . In 
path I of the IFM two RF spin flippers were placed in a 
row, the first one RF 2 " operating at double frequency 2w 
and performing a tt flip, the second one RF"^ 2 operating 
at the frequency u> and performing a tt/2 flip, whereas in 
path II one RF flipper RF£ was placed, operating at the 
frequency u> and performing a tt flip (see Fig. [TJ. RF 2 " 
was equipped with its own Helmholtz coil pair, which en- 
abled it to operate at double frequency. Due to a new 
coil cooling system, thermal disturbances from the coils 
in the IFM were negligible. 




FIG. 1: Experimental setup (see details in text). 



For the W- and the GHZ-like state we entangled 
the path, spin, and energy degrees of freedom of sin- 
gle neutrons. The state was generated as follows: 
The spin-up polarized neutron beam |f 0) entered the 
IFM and split into two partial beams -j= |I t 0} and 

^- | II t 0). In the ideal case, RF 2 " changes i |I t 0) 

to ^75 el7r ^ 2 1 I —2hcj), which then the subsequent RF^ 2 

changes to ^je 47r / 2 |I | -2hu) + ^e i7r |I | -fiw)> while 

RF£ changes 4= |II t 0) to -j=e 17r / 2 |II I -Hu). For the 



symmetric (i.e., balanced) W-state (\W sym )), with the- 
oretical nominal amplitudes a = b = c = -^=, a 
50% stochastic absorber (an indium plate with thick- 
ness 0.75mm) was placed in path II, whereas without the 
absorber an asymmetric (i.e., unbalanced) W-like state 
(I Wasym}) w &s generated with theoretical nominal ampli- 
tudes a = A= and b = c = 4= . For the preparation of the 
GHZ-like state we expect theoretical nominal amplitudes 
d = c = -j=. Here, only RF£! was switched on for state 

preparation while RF 2 " and RF£! , 2 did not operate. 
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Measurement 

Only the O-beam was used for the measurements, 
where a spin analyzing super-mirror (transmitting up- 
spin neutrons only) together with a Larmor accelerator, 
an RF flipper RF" nd , operated at the frequency uj, and 
a DC flipper DC % /2, performing a ir/2 flip, enabled the 
selection of neutrons in dependence of their polarization 
and their energy for detection. Neutrons with the se- 
lected spin properties were counted in the subsequent O- 
dctcctor. 

A full tomography of the generated states would have 
been too complex and partly even infeasible (because cer- 
tain coherences cannot be achieved). As entanglement 
witnesses do not require a full state tomography they 
clearly provide an advantage in this case. 

For the W-state, Eq. (Tl]), the terms (101|p|101), 
(011|p|011), (002|p|002), KOlllpllOl)!, 1(002^1101)1, and 
(002^1011)1 were determined. The terms (101|p|101), 
(011|p|011), and (002|p|002) were determined by inten- 
sity measurements while blocking one path in the IFM 



(the Larmor accelerator and DC„./2 did not operate). 
For the measurement of (101|p|101), path I of the IFM 
was blocked while in path II the beam was manipu- 
lated by RF", and the analysis setup only consisted of 
the super-mirror. The reference measurement was ob- 
tained while RF^ was switched off. For the measurement 
of (011|p|011) and (002^1002), path II of the IFM was 
blocked while in path I the beam was manipulated by 
RF^ and RF^y 2 . Spin- up, spin-down and a reference 
measurement were made. For the spin-down measure- 
ment, RFg nd performed a tt flip, while at the two other 
measurements it was switched off and for the reference 
measurement also RF;;" and RF",„ were switched off. 
The mathematical relations between the measured inten- 
sities and (101|p|101), <011|/o|011>, and (002|p|002) can be 
derived from (fTJ by considering the manipulations of the 
neutron beam in the analysis setup by the beam blocker, 
RF" nd , and the super-mirror (here we had to assume that 
all devices manipulating the beam line were operating as 
previously determined). Each measurement was carried 
out 10 times in order to reduce statistical errors. 




FIG. 2: Typical sinusoidal oscillations of the count rates: a) and b) show contrast scans for |(011|p|101)| and |(002|p|101) 
respectively, together with the reference oscillations (dashed curves) and c) illustrates a contrast scan for |(002|p|011)|. 



In the experimental setups for the determination of the 
cross terms |(011|/?|101)|, |(002|p|101}|, and |(002|/9|011)| 
both paths of the IFM were open and we determined 
these terms by means of the contrasts C of three sinu- 
soidal oscillations while varying the relative phase be- 
tween coherent terms. Coherence between each two of 
the three terms of ([T]) was achieved by three different 
analysis setups. The determination of |(011|p|101)| and 
(002|p| 101) | requires the coherence of terms that are gen- 
erated in different paths of the IFM, therefore contrast 



scans were made by successively turning the phase shifter 
to successively change the relative phase \ between the 
respective two terms. The two terms for |(011|p|101)| 
have the same energy but different spins, therefore in the 
analysis setup DC^ was switched on so that the two 
resulting spin- up terms could interfere 31| . For the mea- 
surement of | (002|p|101) | two terms have to interfere that 
also have an energy difference of Huj, therefore also RF^ nd 
performed a 7r/2 flip, which transformed the two terms in 
question into energetically equal and after T)C„/2 into 
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coherent ones. At each point of these contrast measure- 
ments also a reference point was recorded, where the three 
RF flippers in the IFM were switched off, from which 
reference contrasts C rcf were obtained (the average of 
the reference contrasts amounts to 0.455). The contrasts 
were divided by the corresponding reference contrasts in 
the measurement evaluation. (002|p|011) |, whose com- 
ponents were generated in the same IFM path and had 
different energies, was determined by changing the com- 
ponents' relative phase <f> by successively changing the 
strength of the magnetic field in the Larmor accelerator 
and then making the components coherent by RF" nd per- 
forming a 7r/2 flip. The mathematical relations between 
the measured contrasts and the three cross terms can 
be derived from JTJ) by considering the manipulations of 
the neutron beam in the analysis setup: 2|(011|p|101)| = 
C ab /Clf, 1(0021^101)1 = C ac /CZ f - K0HH101)|, and 
21(0021/31011)1 = Cbc- Each measurement was carried out 
4 times in order to reduce statistical errors and the con- 
trasts and their errors were determined from least squares 
fits. Typical contrast scans are shown in Fig. [5] 

We carried out these measurements without and with 
the 50% stochastic absorber in path II, yielding the asym- 
metric and the symmetric W-like state, respectively. In 
the first case, each point was measured for 60s, whereas 
in the latter case for 80s, to compensate for the total in- 
tensity loss due to the absorber. To test the robustness 
of the generated states, all these measurements were re- 
peated under degrading of coherence between the beams 
in path I and II. The loss of coherence was achieved by a 
shift of the wave packets relative to each other due to an 
aluminium plate placed in path II. For the GHZ-like state, 
©, the terms (101|p|101), (010|p|010), and |(101|p|010)| 
were measured. The determinations of these terms are 
analogous to those described above. 



RESULTS: DETERMINATION OF 
ENTANGLEMENT WITNESS 

The fidelities achieved for the GHZ-like state and the 
symmetric and asymmetric W-like states were 0.985 ± 
0.011, 0.987±0.029, and 0.948±0.022, respectively, which 
values can be mainly attributed to the imperfect efficien- 
cies of the spin flippers. The violation of the inequalities 
and (H]) listed in Tab. U clearly demonstrates the pres- 
ence of genuine multipartite entanglement. The fideli- 
ties achieved for the respective degraded W-states were 
0.646 ± 0.027 and 0.611 ± 0.021. 

Studying the influence of noise, i.e., decoherence, on 
the entanglements is of importance. A wave-packet 
shifter, which is in practice a parallel-sided plate, was 
employed for this study: the shift of the wave pack- 
ets was accurately controlled so that we could control 
the amount of the interference effect in the IFM. The 



State Noiseless With noise 

\GHZ) Ighz = 0.49 ± 0.01 NA 

\W aym ) Iw = 0.47 ±0.03 I w = -0.04 ± 0.02 
\W asym )\l w = 0.46 ± 0.02 \lw = -0.01 ± 0.01 

TABLE I: The almost maximal violation clearly indicates 
multipartite entanglement and thus demonstrates the high 
amount of control in neutron state preparation. In the pres- 
ence of noise (degrading of coherence) no multipartite entan- 
glement is detectable. 



State 


Noiseless 


With noise 


\GHZ) 


h~se V = 0.49 ±0.01 


NA 


\w sym ) 


h- sep = 0.65 ± 0.02 


h-sep = 0.31 ±0.01 


| Vt^ asyrn) 


h-sep = 0.45 ±0.01 


J 3 - sep =0.11 ±0.01 



TABLE II: Despite the strong noise and (potential) loss of gen- 
uine multipartite entanglement the data show that the states 
are not fully separable. This can in this specific case also be 
seen by looking at the fidelity witness. 

situation where all genuine multipartite entanglement is 
lost with retaining bipartite entanglement in the system 
was accomplished and investigated. The analysis was 
achieved via inequalities Ik- S ep[p] < [explicitly writ- 
ten in Eq. ©], which are satisfied for any fc-separable 
state derived in Ref. Q. Thus these inequalities enable 
us to reveal more refined structures. For example, if one 
finds a violation of the inequality Iz-sep[p] for a tripar- 
tite state p, then the state is not fully separable, namely, 
entangled. However, it may still be "only" bipartite en- 
tangled, while a violation of /2-sep shows that the state is 
genuine multipartite entangled (while the criteria Ighz 
and Iw are optimized to detect the different types of gen- 
uine multipartite entanglement). The results are listed in 
Tab. mi The measured values are greater than and thus 
the non full-separability is proven, even in the noisy case. 
The measured values are greater than and thus the non 
full-separability is proven, even in the noisy case. (The 
fidelities also prove inseparability in these systems, albeit 
with a smaller experimental significance). 



In the present neutron interferometer experiment spe- 
cific tripartite-entangled states have been successfully 
generated with high fidelity and their entanglement con- 
tent quantified using nonlinear witnesses. This gives a 
basis for studies on the mixture and transition between 
these states and the totally decohered state, by which 
characteristic properties of tripartite entanglement can be 
demonstrated. In such studies, influences of various kinds 
of decoherence on entanglement will be investigated: in 
particular, a method involving controllable decoherence 
(e.g. like in experiment [32}) will play an important role. 
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CONCLUSIONS 

We have succeeded in generating tripartite entangled 
W-like state families and a GHZ-like state with high fi- 
delity within neutron interferometry. Genuine tripartite 
entanglement of these states is analyzed by proper non- 
linear witnesses. It is demonstrated that, only by measur- 
ing few matrix elements, different types of multipartite 
entanglement can be proven to be present in a single- 
neutron system. Moreover, via the control of the system 
under investigation we could decohere some components 
of the system and thus demonstrate the loss of multipar- 
tite entanglement while inseparability is still kept. These 
results show that quantum correlations can indeed man- 
ifest in different subsystems of a single physical entity. 
Our revealed substructure of an entangled single-neutron 
system and high controllability opens the possibility for 
testing more fundamental and sophisticated properties of 
neutrons via exploring just this substructure, e.g. its re- 
lation to topological phases [33l435| or to the effect of a 
gravitational field 
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